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Abstract. We study thermal emission of a cavity quantum electrodynamic system 
in the ultrastrong-coupling regime where the atom-cavity coupling rate becomes 
comparable the cavity resonance frequency. In this regime, the standard descriptions 
of photodetection and dissipation fail. Following an approach that was recently put 
forward by Ridolfo et al.|arXiv:1206.0944 , we are able to calculate the emission 
of systems with arbitrary strength of light matter interaction, by expressing the 
electric field operator in the cavity-emitter dressed basis. Here we present thermal 
photoluminescence spectra, calculated for given temperatures and for different 
couplings in particular for available circuit QED parameters. 
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1. Introduction 

The quantum theory of photodetection as originally formulated by Glauber pLj is 
a landmark for quantum optics and has occupied a key role in understanding 
radiation-matter interactions. Recently, a new regime of interaction cavity quantum 
electrodynamics (cavity QED) has been reached experimentally where the coupling 
strength between an emitter and the cavity photons becomes comparable to the 
transition frequency of the emitter or the resonance frequency of the cavity mode 
El H la El E]. In this so called ultrastrong coupling (USC) regime [HI El HOl [21 |3l E] 
the usual rotating wave approximation is no longer applicable. As a consequence the 
number of excitations in the cavity-emitter system is no longer conserved, even in the 
absence of drives and dissipation. Hence in this regime, the standard descriptions of 
photoluminescence of a strongly coupled system fails [121 [131 According to the 
very recent result [15], it is possible to calculate the correct output radiation avoiding 
unphysical contributions to the emission from the vacuum due to the presence of the 
counter rotating terms. Here we present the output thermal emission for a cavity QED 
system where a two level system (TLS) interacts with one cavity mode in the ultrastrong 
coupling regime. In particular, we analyze the case in wich the interaction has the form 
of the Rabi Hamiltonian, showing how the spectra modify their shape as the coupling 
strength or the thermal feeding increases. Our results can be measured in a wide recent 
experimental setups [2l|3l [H [Sj [6l [7] . 

2. Input-output relations 

Applying Glauber's idea of photodetection, we here introduce a full quantum theory to 
study the thermal emission in the USC regime [I5]. This requires a proper generalization 
of input-output theory [16], since the standard relations would for example predict an 
output photon flux that is proportional to the average number of cavity photons, i.e. 
('^out(^)'^out(i)) {a){t)a{t)) for vacuum input. Hence an unwary application of this 
standard procedure to the USC regime would predict an unphysical stream of output 
photons for a system in its ground state which contains a finite number of photons 
due to the counter-rotating terms in the Hamiltonian. It has been shown [15] that 
applying the Glauber's formulation of photodetection [Ij, it is possible to derive n-t/i 
order correlation functions for the output fields which are valid for arbitrary degrees 
of light-matter interaction, by expressing the cavity electric-field operator in the atom- 
cavity dressed basis. An ideal detector absorbs a photon with a probability per unit time 
that is proportional to {E~ {t)E~^{t)) where E^{t) are the positive and negative frequency 
components of the electric field operator of the output field [H HT]. In the circuit QED 
the same quantities are measured with output voltages which are proportional to the 
electric fields. Following [13], the input-output relations for a cavity that is coupled to 
a one-dimensional output waveguide via an interaction between the cavity field X and 
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the momentum quadratures 11^ of the waveguide field outside the cavity are 

aout{t) = ain(t) - 2— =^^=X+, (1) 

where ec is a couphng parameter and eo is a parameter describing the dielectric properties 
of the output waveguide, v is the phase velocity. The input (output) field operators 
flout(in)(^) are defined as 

aout(in)(t) = ^ / c?wv^e-'"(*-*)a^(t), (2) 
V In J 

where t = ti — )■ +oo for the output field and t = to ~^ ~oo for the input field and 
ai_j{alj) annihilation(creation) operators of the fields outside the cavity. In this way, 
^(^^out(^)'^out(^)) is proportional to the measured {E~ {t)E~^{t)) and describes an energy 
fiux associated to the output light. The standard definition of output fields as used in 
many textbooks, c.f. [18], is recovered if all frequencies of the field are very close to a 
carrier frequency cJ and one may approximate y/co ^ in the integral kernel which 
makes the observed energy fiuxes proportional to photon number fiuxes. One thus 
needs to find the positive frequency component of X according to its actual dynamical 
behavior, c.f. [19], to compute the proper output fields. We do this by expressing X in 
the atom-cavity dressed basis. It is worth to notice that in the USC regime, the positive 
frequency component of X is not proportional to the photon annihilation operator a. 



3. Dynamics of the open quantum system 

We consider the Rabi model, that consists of a linear coupling between a single cavity 
mode and a two level system (TLS). In the following we set for sake of simplicity h = 1 
and also the Boltzmann constant ks = 1- The Rabi Hamiltonian reads 

Hs = cuoa^a + u^cr^cr^ + g{a + a'')a^ (3) 

where c<;o(x) is the bare energy of the cavity mode (TLS), g is the coupling strength 
and (Jx is the standard Pauli operator. Fig. [l] shows a plot of the spectrum of the 
eigenvalues as function of the coupling g calculated for the parameters uq = = I- 
In particular, it is worth to notice that for a given value of g, the system shows 
the characteristic ladder configuration like in the standard Jaynes-Cummings model. 
The possible transitions between the rungs are given by the selection rules due to the 
Hamiltonian structure, that enables these transitions owing to the non-zero values of 
the elctric field components. In order to describe a realistic system, dissipation induced 
by its coupling to the environment needs to be considered. Then, owing to the very high 
ratio g/ojQ, a standard quantum optical master equations fails as it would for example 
predict that even zero temperature environments could drive the system out of its ground 
state. An adequate description of the system's coupling to its environment requires a 
perturbative expansion in the system-bath coupling strength. In order to perform this 
expansion we write the Hamiltonian in a basis formed by eigenstates \j) of Hq, denote the 
respective energy eigenvalues by huj, i.e. HqU) = huj\j), and derive Redfield equations 
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Figure 1. (color online) Energy spectrum of Hq as function of the coupling strength 
for the Rabi Hamiltonian. This plot is calculated for = wq. The level structure 
is analogous to that of the JC model. The ladders correspond to a coupling strength 
g = 0.2 wq. The arrows indicate possible transitions of radiative decay due to the 
Hamiltonian structure. 

[20] to describe the dissipative processes [21] • We choose the labehng of the states |j) 
such that Uk > ojj for k > j and focus on a single-mode cavity with T 7^ temperature 
for the environment. Generahzations to a multi-mode cavity are straightforward. For 
our porpouse, we neglect small Lamb shifts and dephasing contributions as they do not 
alter significantly the output photon, allowing a simpler and lighter theoretical setup. 
We thus arrive at the master equation, 

p{t) = i[p{t),Hs] + Cap{t) + C^p{t). (4) 

The expressions Ca and Cx are Liouvillian superoperators describing the 
losses of the system where C,p{t) = E,fc>,rf(l + n,{Akj,T))V[\j){k\]p{t) + 
E,», rf ne(Afc„T)I)[|A;)(j|]p(t) for c = a, a" and V[0]p = \{20 pO^ - pO^ O - O p) 
and T is the temperature of the thermal bath. Here nd^kjiT) is the number of ther- 
mal photons that feed the system acting on all the possible \k) — )■ \j) transitions. 
Standard dissipators are recovered in the limit (7 — 0. The relaxation coefficients 
V^J^ = 2'7idc{Akj)al{Akj)\Cji^\'^ depend on the spectral density of the baths dc{Akj) 
and the system-bath coupling strength adAkj) at the respective transition frequency 
Akj = ojk — ojj as well as on the transition coefficients Cjk = —i{j\{c~ c^)\k) (c = a, a~). 
These relaxation coefficients can be interpreted as the full width at half maximum of 
each \k) — )■ \j) transition. Since we consider a cavity that couples to the momentum 
quadratures of fields in one-dimensional output waveguides, it is possible to show [15] 
that the relaxation coefficients reduce to F^'^ = 7c ^ IC'jfcP; where 7c are the standard 
damping rates of a weak coupling scenario. 
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4. Results 

According to the input-output relation ([T]) the output spectrum of light is, for input 
fields in vacuum, equal to, 

oc lim 23ft / + r))e*^"dr (5) 

i.e. the Fourier transform of the two time correlation function. In Eq. (|5]) X = 
— iXo(a — a^) and 3ft denotes the real part. In this way, we are calculating the 
photoluminescence of an interacting system that is thermalized at temperature T and 
emits into the vacuum, i.e. into a reservoir at T = temperature. To separate X in 
its positive and negative frequency components, and X~, we expand it in terms of 
the energy eigenstates |j) and find X^ = — X]j fc>j ^fci^jfcb)(^l) where Xj^ = {j\X\k) 
and X~ = {X~^y. Note that X+|0) = 0, for the system ground state |0) in contrast 
to a|0) 7^ 0. For a weak excitation density S{u) is proportional to the emission. Then, 
once the master equation ^ is numerically solved, one can easily calculate the thermal 
spectrum ^ applying the quantum regression theorem [18] . Fig. [2] shows a plot of 
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Figure 2. (color online) Thermal emission spectra for different temperatures: T = 
O.lwo (blue-dotted line), T ~ 0.15wo (green-dashed line), T = 0.2ujo (red-continuous 
line). The other parameters are uj^ = wq, 7a = 7a; = 5 x 10^^ luq, and g — 0.1 loq for 
all three cases. 



thermal emission spectra for different temperatures. The used parameters are scaled 
by the bare mode energy of the cavity and in this figure g = 0.1 uq. If we choose 
uq = 2tt X 10 GHz, the parameters belong to the common range of values of circuit 
QED systems and the corresponding temperatures fall in the region of hundreds of mK. 
In fact, for = 2tt x 10 GHz, the respective temperature for the blue-dotted curve 
is 48 mK, for the green-dashed 72 mK, and for the red-continuous curve is 96 mK. In 
Fig. [3} instead, we increased the value of the coupling up to g = 0.2 uq, leaving the 
other parameters the same as in Fig. [2j One can see that the mean peaks, that in Fig. 
[2] correspond to the values uj/loq ~ 0.9 and uj/uo ~ 1.1, are now separed roughly by 
twice the energy with respect to the previous case. This happens because the energy 
separation of the first two eigenmodes of (Isl), also in this range of interaction, growths 
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Figure 3. (color online) Thermal emission spectra for different temperatures: T = 
O.lwo (blue-dotted line), T = 0.15wo (green-dashed line), T = 0.2u!o (red-continuous 
line). The other parameters are the same used in Fig. [2] except for the g coupling that 
here is g = 0.2 wq. 

linearly with respect to the coupling strength. Looking at the figures |2] and [3] we can 
see that if the temperature increases, the number of resonances present in the system 
increases as well. Indeed, whenever the temperature rises, this increases the thermal 
occupancy nc{Akj,T) allowing to exicite the \k) — )■ \j) transitions. A characterization 
of the resonances, e.g. in Fig. [3} is really straightforward if one looks at the ladder 
scheme (see Fig[T]). Sorting the energies from the lower value to the higher, we identify 
|4) ^ |3), |2) ^ |1), |4) ^ |2), |5) ^ |3), |3) ^ |1) and |5) ^ |2) as shown by the rows 
in Fig. [T] It is worth to notice that in all the spectra that we presented, there is an 
evident asymmetry also in the absence of detunig, in contrast with the results achievable 
within the standard model for the dissipation and photodetection that naturally cannot 
apply in this regime of interaction. This is due to the fact that in the USC regime, i) the 
effect of the thermal feeding of the reservoir acts differetly on the different transitions 
and, ii) each resonance is characterized by a different damping rate. 

5. Conclusions 

In conclusion, we presented a full description of the thermal emission in the USC 
regime, valid for arbitrary light-matter couplings. We showed thermal emission spectra 
calculated for available circuit QED parameters. These results show that the recently 
proposed correlation function are also able to describe correctly incoherent light emission 
from USC systems avoiding unphysical emission from the ground state. Generalizations 
of our study to multicavity devices p2| and to three-level in the USC regime [23112^125] 
would form interesting perspectives for future research. 
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